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An elongation method based on ab initio quantum chemistry approaches is presented. It allows to 
study electronic structures and coherent electron transportation properties of single-walled carbon 
nanotubes (SWCNTs) up to 22nm in length using the hybrid density functional theory The 22nm 
long SWCNT, consisting of more than ten thousands electrons, is the largest carbon nanotube that 
has ever been studied at such a sophisticated all-electron level. Interesting oscillating behaviour of 
the energy gap with respect to the length of the nanotube is revealed. The calculated current- voltage 
characteristics of SWCNTs are in excellent agreement with recent experimental results. It confirms 
the experimental observation that a 15nm long SWCNT is still largely a ballistic transport device. 
The proposed elongation method opens up a new door for the first principle study of nano- and 
bio-electronics. 
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Carbon nanotubes (CNTs) are probably the most stud- 
ied nanomaterials in the last decade, owing to their great 
physical and chemical properties. Among many excit- 
ing applications, CNTs are believed to play an impor- 
tant role in the future electronics [J S @, H 0. This 
view is further enforced by an exciting recent develop- 
ment, namely the utilization of 10 to 50 nm long single- 
walled carbon nanotubes (SWCNTs) 0, 0. It has been 
demonstrated that the electron transport in SWNTs with 
finite length less than 15nm immunes from the opti- 
cal phonon scattering and thereby exhibits nearly ballis- 
tic behaviour at high biasesjj, P However, the finite- 
length SWNTs present a great challenge for the first 
principle theoretical modelling because of the involve- 
ment of vast number of electrons. Over the years, dif- 
ferent approaches, such as semi-empirical Ifj [7j, ab ini- 
tio tight- binding[§ H GU [U [H, G]| , and "integrated" 
[lil lldi Hrl H?| methods, have been developed to describe 
the nano-sized systems, but all suffer from a relative low 
level of accuracy. For such large systems, the most accu- 
rate yet feasible approach is probably the density func- 
tional theory (DFT). In this letter, we present a straight- 
forward elongation method in conjunction with the mod- 
ern quantum chemical density functional theory calcula- 
tions that allows to effectively treat very large nano-scale 
periodical systems without losing the accuracy. This 
method has applied for studying the electronic structures 
and the coherent electron transportation properties of 
SWCNTs of different lengths. The largest SWCNT re- 
ported here is about 22 nm long in length which consists 
of ten thousands of electrons and is described by fifteen 
thousands gaussian basis functions. 

Our elongation method is based on a simple fact that 
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for a large enough finite periodic system, the interaction 
between different units in the middle of the system should 
be converged, and consequently those units in the middle 
become identical. It is thus possible to elongate the ini- 
tial system by adding the identical units in the middle of 
the system continuously. This can be easily done when 
the Hamiltonian of the system is describe in the site- 
representation. Obviously, the precondition for using the 
elongation method is to obtain an initial Hamiltonian in 
the site-representation with identical middle parts. For- 
tunately this condition can be achieved routinely with 
the modern quantum chemistry programs. The elon- 
gation method should be as accurate as the quantum 
chemistry method used for the initial system. We will 
use the (5,5) metallic SWCNTs to demonstrate the per- 
formance of the elongation method as implemented in 
the QCME code[I3. With the same code, the coherent 
electron transport in the SWCNTs is also calculated us- 
ing the generalized quantum chemical Green's function 
method HlGil. 

We first calculate two (5,5) metallic SWCNTs with 
190 (19 units,CNT19) and 310 (31 units, CNT31) car- 
bons, respectively, at the hybrid density functional the- 
ory B3LYP level with STO-6G basis set using GAUS- 
SIAN03 program 20]. The both ends of the finite-sized 
tubes are terminated by hydrogen atoms. The calcu- 
lated Hamiltonian of the CNT19 is used as the base 
for the elongation method to construct a Hamiltonian 
of the CNT31. The calculated highest occupied molec- 
ular orbital (HOMO) and lowest unoccupied molecular 
orbital (LUMO) energies from the elongation method are 
found to be -2.000 and -1.538 eV, respectively, in excel- 
lent agreement with the values of -2.070 and -1.598 eV, 
respectively, obtained from GAUSSIAN03 program. The 
derivation for the LUMO-HOMO gap E g from these two 
different approaches is as small as 0.009 eV. The calcu- 
lated transition probability spectra above the Fermi level 
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FIG. 1: (A) Transition probability (on a log scale) spec- 
trum above the Fermi level, and (B) I-V characteristics of 
the CNT31 (310 carbons, 3.7nm in length) connected to two 
gold electrodes obtained from Gaussian03 (solid lines) and 
elongation (dashed lines) methods. 



for the CNT31 connected to two gold electrodes, as well 
as the corresponding current- volt age curves, from both 
methods, agree very well as demonstrated in Figure 
It should be mentioned that the coupling coefficients be- 
tween electrodes and CNTs is approximated by a sim- 
ple Au-C diatomic interaction. The coupling between 
electrodes and CNTs is a complicated issue and will be 
discussed elsewhere. 

We are thus ready to explore the elongation method 
for even longer CNTs. To maintaining higher accuracy, 
we have used CNT21 (210 carbons, 2.6nm long) tube 
calculated at B3LYP/6-31G level with GAUSSIAN03 
to construct a series of SWCNT systems with units of 
n = 9 + 12 x i (i=2, 3,..., 14) whose length goes from 
4.1 nm to 21.8 nm. In Figure El (A), the molecular or- 
bitals (MOs) around the energy gap are displayed. As ex- 
pected the density of states (DOS) gets higher for longer 
SWCNT. The most interesting observation is that the 
energy gap (E g ) oscillates periodically with the increase 
of SWCNT length. For short SWCNTs, it is known that 
the Eg should oscillate with a period of 3 units jH [H . 
This behaviour is nicely reproduced by our calculations, 
see Figure EtC). The trend given by the short SWCNTs 
seems to indicate that the gap should continuously de- 
crease with the increase of the SWCNT length. However, 
our calculations present a very different picture. For in- 
stance, the energy gap of the CNT33 (0.598 eV) is al- 
ready larger than that of the CNT21 (0.492 eV). This 
observation is further confirmed by Gaussian03 calcula- 
tions with STO-6G basis set. The evolution of energy gap 
with respect to the tube length is clearly shown in Fig- 
ure E[B). It is noticed that the width of the oscillation 
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FIG. 2: Molecular orbit als of SWCNTs with N = 9 + 12 x i 
(i=0, 1, ... ,14) units, corresponding to 1.1 nm to 21.8 nm 
in length (A). The results for the CNT9 and the CNT21 are 
calculated with GAUSSION03 while others are obtained with 
the elongation method. Energy gap E g as a function of the 
units (N) is plotted for long-nano-scale (B) and short-nano- 
scale (C) SWNTs, respectively. 



peak covers 36 units from the CNT81 to the CNT117, 
while the one after it is about 60 units from the CNT117 
to the CNT177. It seems to imply that the oscillation 
of the energy gap with respect to the length of the tube 
will be slowly faded out. Furthermore, discrete molecu- 
lar orbital distribution is found for all the tubes under 
investigation, indicating that these tubes can be quan- 
tized electronic devices. One can at least conclude that 
a (5,5) SWCNT with length of 22 nm is still far from a 
real bulk material. 

By using the QCME code, we are able to calculate 
the I-V characteristics of the tubes with different lengths. 
The electrode structure shown in Figure[2{A) is the same 
as that reported in the experimental work of Javey et 
al0. The I-V curves of short tubes up to the CNT45 
(5.5nm) show stair-like profile in which each step repre- 
sents an opening of a new conducting channel. These 
small devices are all semi-conductor-like, i.e. the elec- 
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FIG. 3: (A): Schematic draw of a SWCNT device with gold 
electrodes; (B) Calculated current-voltage characteristics for 
the CNT9 (l.lnm), the CNT21 (2.6nm), the CNT33 (4.1nm), 
and the CNT45 (5.5nm); (C) Calculated current- volt age char- 
acteristics for the CNT69 (8.5nm), the CNT93 (11.4nm), the 
CNT117 (14.4nm), and the CNT177 (21.8nm); and (D) Ex- 
perimental I-V curve for a 15nm SWCNT0. 



tron stars to flow only after adding certain external bias. 
When the length of the tube gets longer, the DOS of 
the system becomes denser, the steps in the I-V curve 
smear out, and the device shows a metallic behaviour. 
We have compared our calculations with the experimen- 
tal result y . As indicated in Figure O^D), our calculated 
I-V curve for a 14.4nm tube (CNT117) agrees surpris- 
ingly well with the experimental result of a 15nm long 
tube for both the shape and the magnitude. It thus con- 
firms the experimental observation that a device with a 
15nm long tube behaves like the ballistic conductor 4]. 

It might be anticipated that the oscillation of the en- 
ergy gap could result in a similar behaviour for the elec- 
tron conductivity. Indeed, such behaviour is clearly ob- 
served for the conductance of different SWCNTs, see Fig- 
ure 0] The period for the oscillation is also the same as 
that for the energy gap, i.e. 3 units. However, the phase 
of these two oscillations is slightly shifted. The tubes 
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FIG. 4: The conductance for SWCNTs with units N =8, 9, 
21 biased at their first conduction state. 



with 3n + 1 units have the lowest conductance, while the 
tubes with 3n units give the smallest energy gap. 

Another interesting quantum effect related to the 
SWCNT is the periodical distribution of its molecular 
orbitals along the tube axis. Venema et al. 24] car- 
ried out the position-scanning I-V measurements for a 
finite-length metallic SWCNT. The experimental setup 
is schematically recaptured in Figure together with the 
experimentally measured result. In this experiment, a 
STM tip scanned through the top of a 30nm metallic 
SWCNT deposited on a gold surface to measure the cur- 
rent at each positions. It was found that the conductance 
of this tube shows a periodic behaviour with a period 
of about 0.4 nm, roughly covering 3 layers. With the 
elongation method, we have simulated the STM experi- 
ment for two tubes, the CNT33(4.1nm) and the CNT117 
(14.4nm) at external biases of 0.3V and 0.1 V respec- 
tively. The results are also shown in Figure (SJ As 
one can see, the oscillating behaviour is well reproduced, 
which shows a period of 3 layers. It should be stressed 
that our calculations describe only the coherent electron 
transport, while for a 30nm long SWCNT phonon cou- 
pling can become important which might explain why 
the periodicity of the experimental result is less perfect. 
Furthermore, there is a very interesting wave-packet-like 
behaviour in the position depended conductance of the 
CNT117. It seems to suggest that an electron standing- 
wave can emerge in SWCNT longer than 15 nm. 

In summary, we have proposed an elongation method 
that is capable of calculating electronic structures and 
electron transportation properties of periodic nano- sized 
systems at the hybrid density functional theory level. 
With this method, we have studied the length depen- 
dence of the energy gap and the conductivity of the finite- 
length single- walled carbon nanotubes (SWCNTs) rang- 
ing from lnm to 22nm. The calculated results agree 
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FIG. 5: Schematic draw of STM scanning experiment (A) and 
the corresponding conductance measurement B) for a 30nm 
SWCNT|24j. Calculated position dependent conductance for 
the CNT33 (4.1nm) at external bias 0.3V and for the CNT117 
(14.4nm) at 0.1V. 



very well with the full quantum chemical calculations 
and exiting experimental results. A new interesting os- 
cillating behaviour for the energy gap of the finite-length 
SWCNTs is revealed. It is worth to mention that our 
elongation has also successfully applied to 50nm long 
poly(para-phenylene ethynylene)s and 60 pairs poly(G)- 
poly(C) DNA. The calculated I-V characteristics of these 
systems are also in very good agreement with the corre- 
sponding experiments. The details of these studies will 
be published elsewhere. 
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